Distribution of conduction channels in nanoscale contacts: evolution towards the diffusive limit 
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We present an experimental determination of the conduction channel distribution in lead nanoscale con- 
tacts with total conductances ranging from 1 to 15 Go, where Go = 2e 2 /h. It is found that even for con- 
tacts having a cross section much smaller than the mean free path the distribution tends to be remark- 
ably close to the universal diffusive limit. With the help of theoretical calculations we show that this be- 
havior can be associated with the specific band structure of lead which produces a significant contribu- 
tion of partially open channels even in the absence of atomic disorder. Published in Europhysics Letters, 
http://www.edpsciences.org/articles/epl/abs/2005/ll/epl8748/epl8748.html. 
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Electron transport properties in a quantum coherent con- 
ductor are fully characterized by the set of transmission coef- 
ficients {t„} corresponding to the conductance eigenchannels 
of the system [1, 2]. Clearly, the whole set {%,} cannot be 
extracted from the total conductance which only yields infor- 
mation on %i- The possibility of an accurate determination 
of the individual {t/„} was demonstrated by Scheer et al. [3] 
for the case of one-atom contacts, by analyzing electron trans- 
port in the superconducting state [4]. This technique allowed 
to prove that the number of channels for one-atom contacts 
is basically determined by the valence orbital structure while 
the particular values in the set {t/„} depend also on the contact 
geometry at the atomic scale [5]. 

A natural question that arises is how the set {t„} evolves 
as the size of the contact is increased. One would ex- 
pect that for sufficiently large contacts (where the size is 
larger than the mean free path) the transmission coefficients 
would be distributed according to P(T) = 8 (T — T„)) = 
(G) / [2GqT\/1 — r), i.e. the universal distribution function 
predicted for a quantum conductor in the diffusive regime [6, 
]. Although indirect evidence of this distribution has been 
obtained through shot-noise measurements [8] a direct deter- 
mination of P(T) in this regime is still an open experimental 
challenge. 

In the present work we combine experimental and theoret- 
ical efforts in order to analyze the evolution of the channel 
distribution in nanocontacts as the size of the contact is in- 
creased from the atomic-size limit. We obtain the set {t/„} 
from transport measurements in lead nanocontacts with con- 
ductances ranging from 1 to 15 Go- These results are com- 
pared with model calculations in which the effects of geome- 
try, atomic disorder, and band structure can be included. We 
find that even for small contacts (G ~3-4 Go) the channels 
distribution for Pb is unexpectedly close to the diffusive limit. 
We show that this behavior can be attributed to the particular 
band structure of Pb. 



Highly stable atomic scale contacts are formed using a me- 
chanically controlled break-junction (MCB J) [9]. A notched 
lead wire 99.99% pure is glued on top of a flexible substrate. 
The rupture of the wire is carried out at 1 .6 K in He exchange 
gas, thus preventing the freshly exposed surfaces from oxi- 
dizing. The contact can be reestablished and broken with pi- 
cometer displacement resolution and stability. 

It was shown that in one-atom contacts the {%,} can be 
obtained from the experimental current-voltage characteris- 
tic curve (/V) by fitting to the sum of N one-channel TV 
curves, I(V) = Ym=i i(?n,V) [3, 10], which have been calcu- 
lated elsewhere for arbitrary transmission T [11]. 

Scheer et al. [3, 10] determined the transmission of a small 
number of channels for one-atom contacts using a steepest de- 
scent based method [ I 2]. For larger contacts we have found 
that a simulated annealing algorithm is much more efficient. 
It is faster, a fact that becomes increasingly important as the 
number of channels to fit grows, and has the advantage of 
avoiding getting trapped at local minima. We use as accep- 
tance criterion the root mean square deviation from the mea- 
sured IV curve, X ■ The fit starts from a random initial posi- 
tion in the TV-dimensional space spanned by { T„ } . Then a ran- 
dom walk is performed. At each step the new fit is accepted 
if x 2 is smaller than the previous one. If x 2 is larger the step 
can still be accepted, according to a Boltzmann-type factor, 
with a pseudo-temperature which is lowered in a smooth way. 
/V is chosen large enough for the fit to result in several closed 
channels (t < 0.01). 

We show in fig. 1 a) IV curves for different contact realiza- 
tions. These IV curves exhibit subharmonic gap structure for 
voltages below 2 A je due to multiple Andreev reflections [11]. 
At low voltages the IV curves show a contribution from the 
supercurrent peak and is not taken into account. 

In fig. 1 b) we show the channel transmissions obtained 
from the fit to the IV curves in fig. 1 a). Channel in- 
dexes are ordered in decreasing transmission. An estimate of 
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FIG. 1: (Color on-line). Top panel shows experimental IV curves (o) 
in the superconducting state, for contacts of G = 2,4,6, 10, 15 Go, 
and the fits (lines). The value of the superconducting gap A = 1.35 
meV used in the fit is obtained from an IV curve in the tunneling 
regime. The inset in a) shows the theoretical IV curves for a sin- 
gle channel used in the fitting procedure (taken from [11]). Bottom 
panel: corresponding sets {t„} obtained from the fits. The inset in b) 
shows a schematic representation of the MCBJ technique. 

^norm 

(T) = P(T)Go/ (G) can be obtained from a finite num- 
ber of experimental realizations in the following way. A his- 
togram of the {t,,} values is constructed dividing the trans- 
mission axis into 14 equally sized bins. The T values falling 
within the range of each bin are collected from 20 contact re- 
alizations of a given conductance. The result is normalized 
dividing by the conductance times the number of realizations. 
We find that P n0 vm(T) is bimodal, with a significant number of 
channels with high transmission. 

The channel distribution can be also conveniently presented 
using the well-behaved integrated quantity [13]: 

Q(T) ^ [ T dT' P(T')T' , (1) 
(G) Jo 

which is a regular function of T. Q(T) is the relative contri- 
bution to the conductance of all channels with T < T . 

We obtain an estimate of Q(T) from 20 experimental real- 
izations making a histogram with a bin size equal to 0.01 . The 
T values that are lower than T from all realizations of a given 




FIG. 2: (Color on-line). Estimate of Q(T) for contacts of G = 
2,4,6, 10, 15 Go (black, red, green, blue and cyan respectively). The 
Q(T) corresponding to the diffusive case is shown as a dashed line. 
Inset: estimated P norm (T). 

conductance are collected and added. Then the result is nor- 
malized dividing by £ T„ . For a perfect ballistic conductor this 
function is zero except for Q(\) = 1. The distribution shown 
in fig. 2 is not very sensitive to the size of the contact and 
remarkably similar to the universal diffusive limit. This simi- 
larity is rather unexpected due to the small size of the contact 
with less than 20 atoms in cross section. 

Previous experience for the case of one-atom contacts 
has demonstrated that band structure effects play a major 
role in the determination of the set {t„} for a given ele- 
ment [3, 14, 15, 16, 17]. It is thus interesting to analyze 
whether this is still the case for contacts in the intermediate 
range (larger than one atom but smaller than the mean free 
path). For this purpose we have performed model calculations 
of the set {t/„} for idealized geometries of the nanocontacts. 

As in previous studies we use a self-consistent parametrized 
tight-binding (TB) model in which the main features of the 
bulk band-structure are included [14]. We shall consider 
model geometries like the ones depicted in the insets of fig. 3 
in which the neck of the contact is represented by several 
atomic layers of different cross section grown along the (111) 
direction on a fee lattice. This central region is connected 
to the left and right electrodes represented by perfect semi- 
infinite crystals. The corresponding TB Hamiltonian can be 
written as H = L//u.„ % ;j<V. (/x ,<\./i.o = H c + V L + Vr + H L + 
Hr, where Hi r and H c describe the uncoupled left and right 
electrodes and the central part of the contact respectively; V~l,r 
being a coupling term between the central region and the elec- 
trodes. The matrix elements hij. a ,p, where i,j and a,/3 de- 
sign sites and orbitals respectively, are taken from fits to the 
bulk ab-initio band structure. As a self-consistency condition 
we impose local charge neutrality on each site. Disorder in 
the atomic positions at the central region can be included by 
rescaling the hopping elements according to the distortion of 
the different bond length with respect to the bulk values [18]. 
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FIG. 3: (Color on-line). Distribution of transmission coefficients cal- 
culated for idealized geometries in which the number of atoms in the 
neckregion layers are: a) 18-12-18, b) 18-12-12-18 and c) 31-27-27- 
3 1 . In all figures the full and open circles correspond to Pb and Au 
respectively. In panel d) we show the distribution Q(T) correspond- 
ing to the case of panel c). 



Once the TB Hamiltonian has been built the conduc- 
tance G(E) = GoTr \p (E)f(E)~\ is calculated in terms of the 
matrix elements of the Green function operator G r {E) = 



lim^o [E + ir\ - ft] using : 



m 



2f 1 ] /2 (E)G[ N (E)f^ 2 (E), 



(2) 



where fi » are the matrix tunneling rates connecting the cen- 
tral region to the leads [14]. The transmission matrix Pi can 
then be diagonalized in order to obtain the transmission eigen- 
channels and eigenvalues T„. 

In fig. 3 we show typical results for the transmission sets 
that are obtained for Pb contacts of intermediate cross sec- 
tion. The effect of an specific contact geometry is illustrated 
in the upper panels (fig. 3 a) and b)). As can be observed, 
a characteristic feature of Pb is the appearance of steps corre- 
sponding to partially open channels with similar transmission. 
For comparison we also show in this figure the corresponding 
results for Au contacts in order to remark the differences with 
a monovalent metal. Although the detailed structure of these 
curves is strongly dependent on the specific geometry, the be- 
havior for Pb is clearly in contrast to the one observed in Au 
where partially open channels are rarely observed. These dif- 
ferences can be directly related to the particular band structure 
of both metals. While in the case of Au the channels can be 
associated with a single band of s-character at the Fermi en- 
ergy, in Pb the channels arise from the contribution of several 
bands with sp-character. In fact, for the same geometry the to- 
tal number of channels with significant transmission is larger 
for Pb than for Au. When the cross section of the contact is 
increased (fig. 3 c)) the partially open channels in Pb tend to 
define a rather continuous distribution which fairly resembles 
that of a diffusive conductor. In panel d) we compare the cor- 



: Pb 


0=0 2 ^stti 


- 


. A 1 J/ 1 

a=0.1 /ja 




\ o=0.0 

O U.UJ 
l.l. 


- Au 




■ 


/J f 

*i n 

s J 
si JA 

s J fl\ 

S J / 1 
s J J j 1 

i — 1 




i.i. 



1.0 
0.8 

0)0.4 
0.2 

1.0 
0.8 

0^0.4 
0.2 
0.0 

0.0 0.2 0.4 0.6 0.8 1.0 

T 

FIG. 4: (Color on-line). Effect of atomic disorder in the calculated 
distributions Q(T) for Pb and Au. The curves correspond to different 
values of the disorder parameter a described in text: (black line), 
0.05 (red), 0.1 (green) and 0.2 (blue). The symbols correspond to an 
average of the experimental results from 10 to 15 Go. The dashed line 
corresponds to the universal diffusive distribution. 



responding distributions Q(T) with the result expected for a 
diffusive wire. 

As a result of these previous calculations one can conclude 
that, even in the absence of atomic disorder, Pb contacts de- 
viate strongly from the ballistic behavior observed in the case 
of Au. 

In a second step we have analyzed the influence of atomic 
disorder in the calculated distributions. Disorder is assumed 
to be localized in the neck region (we will discuss later the 
effect of disorder on the bulk electrodes). The degree of disor- 
der is characterized by the mean square deviation of the bond 
lengths with respect to the idealized structure, O". fig. 4 shows 
the evolution of the distribution function Q(T) as a is in- 
creased for a particular contact geometry. As can be observed 
for Pb the distribution tends rather quickly to the diffusive one 
with increasing disorder. We also show the experimental dis- 
tribution in which several cases with total conductances rang- 
ing from 10 to 15 Go have been averaged. This distribution 
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is in fairly good agreement with the calculated one for the 
case a = 0.1, which roughly corresponds to a mean free path 
l e ~ 10 nm. This rather low degree of disorder is consistent 
with what is known both from molecular-dynamics simula- 
tions and images made by transmission electron microscopy 
of metallic nanocontacts [5]. 

It should be stressed that these distributions still exhibit 
some deviations with respect to the universal curve, specially 
for high transmission. We can associate these deviations with 
the fact that the size of the contact is still considerably smaller 
than the mean free path. 

In contrast to the case of Pb, Au exhibits a clear deviation 
from the universal distribution even for the larger values of 
a considered. Notice that a = 0.2 already corresponds to an 
unexpectedly high degree of disorder [19]. 

As a final issue we estimate the effect of disorder on the 
bulk electrodes by means of a phenomenological model. In 
this model the electrodes are diffusive conductors character- 
ized by a classical probability P c \(t) for returning to the con- 
tact region after a given time t. Modelling the electrode as 
a cone with openning angle y this probabilitiy is given by 

Pd(t) = VF/h.\/3nkl (Df) 3/2 (1 - cosy) j, where k F and \> F 
are the Fermi wavevector and velocityrespectively and D = 
vf4/3 is the diffusive constant [20]. Following the arguments 
of ref. [20] the correction to the transmission distribution can 
be written as 

8P(T) = (£8(T- x n ) (1 - T„))p retum , (3) 

n 

where T e = l e /vp is the elastic scattering time and /? re tum = 
J™P c i(f)df. Assuming that typically y~45° and that l e can 
range from ~ 10 nm to ~ 4 nm in Pb, one obtains /7 re tum ~ 
0.005 — 0.02 leading to a rather small correction which hardly 
modifies the theoretical results already shown in fig. 4. 

In conclusion, we have presented a combined experimen- 
tal and theoretical analysis of the conduction channel distri- 
bution of nanoscale contacts with several atoms in cross sec- 
tion. The experimental method demonstrates the applicability 
of the technique developed by Scheer et al. [3, 10] to much 
larger contacts. It is found that for Pb contacts the distribution 
tends rather fast to the universal diffusive limit. We show that 
this behavior can be associated with the specific band structure 
of Pb which favors the appearance of partially open channels 
in contrast to what is predicted for monovalent metals like Au. 
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